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PREFACE 


The present report is part of a continuing effort to find 
improved methods for ground-based detection of nuclear explosions 


conducted at high altitudes and in space. 


SUMMARY 


The electromagnetic signal generated by interaction of 
nuclear—-explosion—induced electric currents with the geomagnetic 
field can serve.as a basis for detecting such explosions at extreme 

| distances from the earth. In the present paper, an earlier theory 
of this phenomenon developed by the authors is extended to allow 
estimates of the maximum detection range. It is concluded that 
this electromagnetic signal may be detectable at distances of the 
order of 10° V Yee kilometers from a nuclear explosion of total 


yield Yeq kilotons. 
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I. INTRODUCTION 


In a previous paper a the authors considered the detection of 
electromagnetic radiation from a nuclear explosion in space caused 
by Compton-recoil electrons from gamma rays scattered in the 
materials of the nuclear device or in a surrounding material shield— 
the so-called case signal. It was concluded that the case signal 


—3/ 2 
had a peak electric field strength of about 10° R 3/ 


vim ata dis— 
tance R kilometers from the explosion. The maximum detectable 
range of this signal was estimated to be i0° km, independent of the 
explosion yield. 

More recently“, the authors considered another mechanism 
for the generation of electromagnetic radiation from nuclear explo- 
sions which might be useful for detection. This mechanism was 
based on the interaction of nuclear explosion—induced electric currents 
with the geomagnetic field. The theory of this mechanism was apphed 
to explosions which take place weil within the earth's atmosphere. 


It was indicated, however, that this same geomagnetic interaction 


mechanism applies to high—altitude explosions where it can lead to 


We. J. Karzas and R. Latter, Phys. Rev. 126, 1919 (1962). 


ow. J. Karzas and R. Latter, Journ. of Geoph. Res. 67, 


4635 (1962). 
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quite intense electromagnetic signals at great distances from the 
burst point-—thereby suggesting its usefulness for detecting nuclear 
explosions in space. In the present paper, this earlier theory of the 
geomagnetic interaction mechanism will be extended so as to permit 
an éstimate of the intensity of these electromagnetic signals observed 
on the earth's surface from nuclear explosions conducted at the great 
distances of interest to the detection problem. It is shown that this 
signal can be detected at significantly greater distances than the case 


signal. 


TI. GENERAL CONSIDERATIONS 


For nuclear explosions below an altitude of 100 or so kilo— 
meters the primary source of the electromagnetic signal is Compton— 
recoil electrons produced by explosion gamma rays interacting with 
the atmosphere. Above 100 or'so kilometers, x rays” from the heated 
materials of the nuclear device generate photoelectrons which also 
contribute a radiating current. 

The currents produced by the gamma rays and x rays from 
high—altitude explosions have maxima in two regions. One region is 
near the explosion where the gamma-ray and x-ray intensities are 
highest even though in this region the air density may be quite low. 
The second region is in the earth's atmosphere at altitudes of about 
20 to 40 or so kilometers for gamma rays and of about 70 to 110 or 
so kilometers for x rays. These altitudes correspond to the region 
of peak absorption for explosion gamma rays which have a mean 
energy of about one mev and for explosion x rays which have a Planck 
distribution corresponding to a nominal temperature of about one kev. 

In our treatment of the detection problem, radiation from the 


currents near the burst point will be neglected. For explosions at 


Sat lower altitudes, the range of the x-ray photoelectrons is 
too short to contribute an appreciable current. Moreover, the x rays 
are confined to a relatively small region about the burst point, which 
region is made conducting by ionization from the earlier emission of 
gamma rays and so cannot radiate. 


altitudes which are not too great, the ionization in the atmosphere 
produced by the gamma rays and x rays is sufficient to absorb this 
radiation before it reaches the earth's surface. For explosions at 
extreme altitudes, where the gamma-ray and x-ray ionization is 


negligible, this radiation can reach the earth's surface, but it may 


still be neglected, since its intensity turns out to be substantially 
less than that from the atmospheric currents. This lower intensity 
results from the very long mean free path for production of currents 
near the burst point, which more than compensates for the closeness 
of this region to the burst compared to the great distance of the 
atmosphere. Moreover, the intensity of the signal from the region 
near the burst point is further degraded by dispersion through having 
to traverse the low density ionized space medium for great distances 
and then having to traverse the high density ionized ionosphere. 
Figure 1 shows schematically the geometry of the explosion. 
Since the currents contributing to the electromagnetic field at the 
observation point are confined to a relatively small region near that 
point, this geometry may be simplified by assuming that the geomag— 
netic field is uniform (See Fig. 2). Spherical polar coordinates will 
be used. The origin will be located at the burst point and the polar 
axis in the direction of the geomagnetic field appropriate to the 


particular current source being considered. 
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ill. ELECTRON CURRENT AND DENSITY 


The number of gamma rays or x rays emitted by a nuclear 


explosion per unit time is 


Nt) == ft), (1) 


where Y is the yield of the explosion in the form of gamma rays or 
xrays, E is the mean gamma-ray or xray energy, and f(t) ex— 
presses the time variation of the gamma-ray or x-ray emission, 


normalized so that 


| dt f(t) =1. (2) 


The rate at which primary electrons (Compton recoils and photoelec— 


trons) are produced at a distance r in direction 6, from the explo— 


sion 1s 
Api Es t)=g(r) f (t-—r/c), (3) 
where 
— {| dr/ 
ay % de (4) 
4aqrr i 


and the integral extends from the burst point to the point under con— 


sideration. \ denotes the mean free path for production of primary 


electrons, 
» FR o /p * 3x 10° p /pcm (5) 
v yo "oO O ; 
h =} / vip? /p cm (6 
x xo Po/P 4x Pf? ? 


where p/p is the ratio of sea—level air density to the air density 
at the point of interest, and E.. is the x-ray energy in kilovolts. 
The primary electrons are slowed down by the atmosphere and are 


stopped in a distance given by their range, R , where 


R = Roy polp ~ 3x 10 p/p cm, (7) 
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~3_2 
R. = Rox poi Pp =~ 10 EL pole em. (8) 


For simplicity we shall assume that the electrons maintain their 
initial speed, vy? throughout their range and then abruptly stop. 

At a particular point in space and time, the primary electron density, 
secondary electron densities, and currents are determined by adding 
the contributions from all primary electrons made earlier and which 
originated at appropriate positions and times. The detailed calcula— 
tion is given in the Appendix. We give here only the results. 


In the absence of a magnetic field, the primary electron 


density at position r and at time ¢t is 


Riv, Vo 
Dri t)* st) | dv'f i T é - 2)| , (9) 
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where t = t-—r/c. If we assume that the secondary electrons are 
made at a uniform rate during the slowing down of the primaries, 


then their rate of production is obtained from the density of primaries 


vq Bs 
- - O = PFl . . a baal 
by multiplying by Rp? where q 3360 and Eri is the initial 
energy of a primary electron. Thus 
o Vot Riv, 
ee t — al — 
Noe (r,t) a et) | dr |. n (1 “9/2 , (10) 


and the secondary electron density is then 


3) Riv, 
n (r,t) ¥ > g(r) | dt | dr''f [n-r'av lel}. (11) 
O 


With no magnetic field, the net electron current is in the radial 
direction and is 
R/ Yo 
j ©; t)+- ev getr)e | art] 7—H-v,/e) | ; (12) 
O 


where ¢ equals the mean cosine of the angle of scatter and is the 


effective fraction of primaries moving in the radial direction. For 


the photoelectric effect, ¢«~ : vie and for Compton scattering 
es]. 

In the presence of a magnetic field, there are two modifica— 
tions to the above results. First, currents are generated in the 


transverse as well as the radial direction. Secondly, as the electrons 


are accelerated by the magnetic field their radial component of 


velocity is reduced which, in the gamma-ray case, causes the 

current to lag in time behind the gamma-ray pulse and which, in the 
x-ray case, causes an additional lag that is a small increment to the 
lag already caused by the low initial photoelectron velocities. The 
modified results, applicable both to the gamma rays and to the x rays, 
are 


Yo T Riv, 2 2. sinwt"' 
tl a i oe i . 
nook t)* 4 R g(r) Jar | dt | (1—Bcos 6)t +8sin 9 , (13) 


Riv, ; 
j ©, t)* —eev g(r) ; att | (1—8cos 6) T' 
O 


1 t 
+ 8 sin”9 sine | cose + sin“@ cOBW “| , (14) 
Riv, 9 
i,{z. t)R eev g(r) | ed ie (1—B cos 6) T' 
oO 


7 | 
+ 8 sine sina sin @cos 8 (coswt'— 1), (15) 


Riv ; 
(r,t) —ecv g(r) | Pant 1 (1-8 cos “9) 1! 
Oo 
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+ Bsin sin @sinwT' , (16) 


2. si 
9 snot] 
i) 
where for gamma rays, ¢« ~ 1 and B = v/e and for x rays, 
€ = = 8, and terms in the brackets involving 8 = viic << ] 

can be neglected. § is the angle between the magnetic field 


| and the direction of motion of the x rays and gamma rays, 


and w is the Larmor frequency for the primary electrons, 


namely, 


» 2 2B = a- ve 4 h¢ a (17) 


IV. ELECTROMAGNETIC FIELDS FROM 


HIGH-ALTITUDE CURRENTS 


Equations (13) to (16) show that a nuclear explosion produces 
a current and electron density in the atmosphere which expands as 
a pulse at the speed of light radially from the source. The magnitude 
of the pulse is proportional to g(r) of Eq. (4) which is large both 
near the source and in the earth's atmosphere, where most of the 
radiation is absorbed. As remarked béfore, we shall limit our 
considerations to the pulse produced in this latter region. The time 
behavior of the pulse is that of a function which varies rapidly with 
the retarded time +t = t—r/ec, where r is distance from the 
explosion. Itis convenient, therefore, to rewrite Maxweli's equa— 
tions in terms of the retarded time. 


In the conventional form, Maxwell's equations are 


3E 
1 ~— Ang 
=o + | 18 
, 3B 
vx ott (79) 


(20) 


V. E = 47p. (21) 


Using the transformations to retarded time T=t—T /c , namely, 
(22) 


— (23) 


where e. is a unit vector directed radially from the explosion point, 


Maxwell's equations become 


aB aE 
1 — 1 — an. 

—_—_ — ees 24. 
VxB- oO S*5n > o BT ve 4? (24) 
1 3B 1 EB 
= 34 VxE+ re xX >>> (25) 

; °B.. 
7. B--— =O, (26) 
oF 
l r 
- Ns 27 
v.E = Ot 4ip (27) 


As usual, j represents the total current, including "back currents" 
or polarization currents in the medium. Conservation of charge, 


which follows from these equations, becomes 


(28) 


It is convenient to eliminate E from Maxwell's equations. 


and get equations involving only B, and conversely . Doing this, we 


obtain the following equations 


2 47 20 ar.) _ 

“WB, Te WXD, Fe a2 ar B,) = 0, (29) 
2 230 [13 Qar . |. 

WB WeDo OF E ar eB) ~ = i, | = 0, 0) 


2 4 230/13 Qn 
-a), ~ Zoxpy+22|2 3 oay+ i] - 0, (31) 


> 


OT or r 
~(V°E), + 4m Wp) + =e [2 Rony + i,|- 0, (33) 
- WE), + Ar Wey + 2 E sr E,) + = i = 0, = (34) 
where 
ae tO gy. (35) 


We now observe that the variation of the currents with dis— 
tance r is much slower than the variation of the currents with cr. 
This latter variation extends over times t much less than a micro— 
second’, corresponding to the emission times of the gamma rays and 
x rays from the explosion. The spatial variation of the currents 


extends over distances like the atmospheric scale height (~ 10 km), 


4 
Latter, R., R. F. Herbst and K.M. Watson, Ann. Rev. of 
Nuc. Sei. II, 371 (1961). 
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which corresponds to times of tens of microseconds. For the transverse 
currents, the resultant electromagnetic fields have a rapidly varying 
time character similar to the current sources’. Thus, the time 
variation of the fields from the transverse currents is much more 

rapid than their spatial variation. We are thus led to an approxima— 

tion in which we keep only the last bracket in Eqs. (30), (31), (33) 


and (34). That is, 


ao11 3a Om. lm 

ST E sr FED * | < oe: (36) 
atl 3a ou 

Oo} fo) + SEs le 3 
OT 3 ar (PE) c i ; (37) 


and similarly for the magnetic fields B. and B Q° For the radial 
equations, to a similar order of approximation, we obtain a non— 


radiating field given by 


1 3 ob. 

—-_ — —— + ~ 

r or|- < OT ni.) 0; (38) 
and 

ta 8B, on a(sin ej ) 

r oar OT r sing a8 


For the radial currents, the radiated field does not have 
a time character similar to the current sources. Rather its char— 
acter is determined by the spatial extent of the current source. 
For the currents generated in the atmosphere, the spatial extent 
of the currents is very great and therefore the radial components 
contribute a very low frequency electromagnetic signal. This signal 
is, in general, weak relative to that due to the transverse currents 
and therefore may be neglected. 
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where we have approximately p* 2 j,. from Eq. (28). Integrating 


these equations, we arrive at the final approximate equations 6 


ob. 

> + 4m}, = 0, (40) 
1 3a 2n 

, OL (rE,) a jg 0, (41) 
1 3 27 

_— + 2S = 

- OF (rE) a Je 0, (42) 


with the corresponding results for the magnetic field. The first 
equation, Eq. (40), is the familiar one for the radial electric field. 
We shall not be concerned with this non-radiating component, but 
only with the propagating transverse fields. 

We recall that the current density j is to include the second— 
ary currents in the plasma as well as the primary electron currents. 
Thus, 


jziv. ti , (43) 


To determine dsec we use the Lorentz model. We assume that the 
plasma electrons are created at rest and that their motion is determined 
only by the local electric field and by their collisions with the air. 

Thus, if x(t) is the displacement of a plasma electron from its initial 
position, which we assume small comparedtoimportant wavelengths, 


Xtv X= E. (44) 


S Results analogous to Eqs. (41) and (42), but in a planar approxi— 
mation where the factors of r are not included, have previously been 
obtained by C. Longmire (unpublished 1964). 
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For an electron created at time t., the current, — ex , is given by 
; ee t _- (t—t') 
— . = —.. t f 
eX(t;t) = Ja tte E(t!) . (45) 
5 


if hoot ) is the rate at which the electrons are made, then the 
current density is given by 


t 
itt) = — [dt AL.) ektst (46) 


or, from Eq. (45), 


| dt'e ° E(t!) no (t!) - (47) 


Clearly, if the plasma electrons were disappearing, as well as being 
created, the same final expression would result— Noo (t) then being 
the net electron density. Thus, we are led to the following result 


for the plasma current 


2 t —v(r) tt-t'] 


M e t f 
plasma (r, y= md dt te E(r,t Neo (st ) . (48) 


If we substitute this expression for the plasma current into 


the field equations, we obtain for E, and E 


p 
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Sr PE ) ° = < far rE,(r, Tne oot, rye. Vv) (tT 1') 
+ =T rj (e,1) =0, 49) 
+ <2 ey E> tT) = 0 (50) 


and similar results for B 9 and B . These equations are not, in 
? 


general, soluble because of the spatial dependence of vy (r) in the 


exponential. 


However, if v,&) is large compared to the frequencies in 


E and Doe? the equations simplify to 


) att ; ae 
sr ft reE +a rj-0, (51) 


where it is understood that E 7 E, and j7j 9° This equation 
? 7D 


is readily soluble and yields 


27 . Ut i" 
— — | ar o(r'", rT) 
co — 


r 
E(r,1) = —- “a 2 | dr'r' j(rine , (52) 
O 
where o(r, t) is given by 
ef 
o (r,t) = n__. (t,t) » (53) 


mV sec — 
C 


In the opposite extreme, for a collisionless plasma with v,@) =0, 
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» @ 


the equations simplify to 


= ane ot ; 
Ar (rE) + | dt'n sect? 1) BET) +> r j(r, rT} =0, (54) 


7 
f which has the solution 
E(r, 1) =— = : fer [ dr' r'! od (r', T') 


1/2 


ore rl T (1 ett 
I, fe je ee (55) 


Analogous results for the magnetic field are readily derived. In 


particular, it is easily shownithat B = Ey and B, =—E§5 
? 2 


Thus, to obtain the transverse electric field from a pulse of 
transverse current moving with light speed, we add all the contri— 
butions along the line from the explosion to the observer. In the 
conducting plasma these contributions are reduced by absorption 


between the source and the observer; in the collisionless plasma, 


dispersion is introduced by the plasma. 


Ue Vv is independent of position, it is easy to show that the 


solution of Eqs. (49) and (50) is obtained from Eq; (55) by replacing 
v v 
E(r,t) with E(r, te “and j(r,t) by j(r, the 
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V. DETECTION CONSIDERATIONS 


It is clear from Eqs. (52) and (55) that the electromagnetic 
field observed from high-altitude nuclear explosions depends only 
upon the distribution of current and electron density along the line 
between the burst point and the observation point. However, because 
of the exponential variation of atmospheric density, the magnitude of 
the electromagnetic field depends quite sensitivity upon just where 
the currents and electron densities are generated in the atmosphere. 
Moreover, the occurrence of the conductivity o(r,T) in the exponen— 
tial of Eq. (52) indicates that the magnitude of the electron density 
at each position in the atmosphere is also quite important. The 
result is that the magnitude of the observed electromagnetic field. 
is quite sensitive to the gamma-ray and x-ray yield of the explosion, 
the distance between the observation point and the burst point, the 
temperature of the explosion, and the angle of incidence of the 
radiation on the atmosphere’. A comprehensive description of the 
electromagnetic field would therefore require a parametric study 
for beyond the intent of this paper which is simply to indicate the 
possible usefulness of this signal for detection. What we shall do 

Stfects due to the earth's ionosphere are not important in 
the present considerations, since the currents which generate the 
electromagnetic signal are produced below the ionosphere. For a 
detailed discussion of the deposition of explosion x rays and gamma 


rays in the atmosphere, see R. Latter and R. E. LeLevier, Journ. 
Geophys. Res. 68, 1643 (1963). 


§ 


here is simply indicate the crude order of magnitude of the electro— 
magnetic field strength under typical explosion conditions and then 
indicate roughly the limiting detection distance. 

First, we shall consider explosions which are not too distant 
from the earth, which have a temperature of about one kev, and 
for which radiation is incident almost normally on the atmosphere. 
Specifically, these conditions mean that the ionization produced 
deep in the atmosphere by the x rays and gamma rays is high, so 
that Vv, in the regions of interest” is large, allowing the approxi— 


mation Eq. (52) to be applied, and, moreover, 


or 
— | dr'o(r,r) > 1. (56) 


oO 


This condition implies® for x rays that a one kiloton explosion 
is within about a few times 10° km from the earth's surface 


2 
and for gamma rays, within about a few times 10 km. 


oor low pressures, large electric fields and electron-volt 


energies, the collision frequency of electrons in air is pproximately 


independent of the E-field and has a value v, = 4x 101 * p/p. sec 
see R. A. Nielsen and N. E. Bradbury, Phys. Rev. 51, 69 (1937). 


To estimate the magnitude of this radiation, we first simplify 
Eq. (52). We note that the factor r' in the integrand of Eq. (52) can be 
replaced by r, the distance to the observation point, since the 
important currents are confined to a thin layer in the atmosphere 
near the observation point. Thus 


r 

r —-— [ ar"ote", T) 

E, (r,1)* -— | dr'j, (r', T)e r! . (57) 
? C 64 > 


Now we note that for r‘ sufficiently near r, 


— [arrote" tT) 0 and j(r',r)*+ 0, (58) 
since the gamma rays and x rays are rapidly absorbed by the atmos— 
phere below some low altitude. As r‘ decreases and the altitude 
increases, we finally reach a radius ry at which the x-ray and 
gamma-ray currents and ionization increase rapidly. Thus we can 
rewrite Eq. (58) as 

r _ on 


r 
; = |°dr"o(e", 7) 


E (r, 7) = _ 2a | Y dr! jy (r', The r! ’ (59) 
8, Cc 5 »® 


or making an obvious change of variables, 
ro 
a PP eta tel) 


a (r(x), 7) ~ 


—X 


E, (r,1)*~ - | dx (60) 
vy 


j (r'(x), 7) 


varies slowly compared to e 
g (r(x), 7) : yfeome 


It is easy to show that 


~2 |— 


and therefore 10 


r 
oO 
5 (ry) — otf arts (e"', 1) 
E (r,t) - "8,97 oO" ‘L-e © 
ap? o(r, 7) C ° ) , 
Jp (r 57) 61) 
o(r_,7T) 


Substituting Eas. (15), (16) and (53) into Eq. (61), we find 


€mw VR Jy (rT) 
E = ——— as a: ; (62) 
8, eq N(t) 
where 
. Riv, 2 2. sinwT' 
N(t) =u | at" | art| —(1-—Bcos 6)T'+8sin 8 ——| > (63) 
— 2 oO 
Riv, 2 2. sinw T' 
J (7) = | ante] (1—f cos §)T'+Bsin 6 sine sinécosé (coswt'—1) , 
) 
IV, 2 2. Sinwt'! 
Jf = drt] a (i1—Bcos 9)T'+6sin §—— sin§ sinwr' . (65) 
<) 


We note that for gamma rays, which deposit most of their energy in 


the altitude range from 20 to 40 km, 


=. ~ 2x 10 to 3x 10° sec. (66) 
Oo 


Equation (61) gives the result obtained by a more direct 
calculation in the paper cited in footnote 2. 


(64) 


ee ee Pee be 


ia. > ‘T-. bie 
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Since the gamma rays are emitted in times of about 10 sec, we see 
that order—of-magnitude—wise f(t) can be replaced by a delta function. 
Similarly, for x rays which are emitted in times* considerably less than 


10° sec, and are absorbed in the altitude range from 70 to 110 km, 


= ~ 10 to 5x 10° sec, (67) 


O 
—2 
where R~ 10 p/ p, cm for electrons with energies of a few kilovolits, 


so that again f(t) can be replaced approximately by a delta function. 


Thus, 
~ wR 
~ oR wR 
in 6 _ 
,(n) ~ — Sin coe 6 (1 208 x) x < -< 
l—Bcos @—Bsin 8cos X 
w KR 
J(n 5 sin cen x . ,x< ee 
? 1-8 cos*@— 8 sin’ @cos X 
7 wR 
= 0 ; a > Bc. ) (70) 
where X is determined from 
2 2, 
(l1—8cos 6)X—fsin OsinX wrt. (71) 


Numerical evaluation of J, pi N> using 8 =v ic = 0.95 correspond— 
3 


ing to 1mev Compton—recoil electrons and Bzv /c =0.05to 0.1 
O 


corresponding to 1 to afew kev x-ray photoelectrons, gives the peak 


yr © 


wk 1.67 x 10° cps, e =4.8x 10. 
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J oi N~ 1/{1-8) . Hence, order—of—magnitude—wise 


emuvR 
Cc 


E ——— 
9, eq(1—) 


(72) 


For gamma rays wherein «x 1, m =9x 10 7° gm, 


10 2 


qa 3x 10°, and R~3x 10 


pi pcm, we have 


E ~ 2esu or 6x 10° vim. (73) 
9,0 | 


2 
And for x rays, wherein?! e~Q0.1,0%5.3x 10° cps, q* 10 , 


—2 
and R* 10 p/p. 


E, ~~ 3.5 x 10° esu or 10v/m. (74) 
7 


Both of these fields are sufficiently intense so as to be easily 


detectable. le 


We now turn to the problem of estimating the intensity of the 


electromagnetic signal from explosions at extreme distances from 


thy mean photoelectron energy of about 3 kev is assumed, 
since the x rays which contribute most to Ey , are those which 
penetrate deepest. For the range of low-energy electrons, see 
Eq. (8) and L. Katz and A. S. Penfold, Rev. Mod. Phys. 24, 28 (1952). 


12s ctually, for the assumed close distance of the explosions 
and the assumed normal incidence of the radiations on the atmosphere , 
the pre—ionization produced by the gamma rays is sufficiently intense 
to absorb the x-ray produced signal before it reaches the earth's 
surface. Thus, under these circumstances, only the gamma-ray 
produced signal is useful for detection. 
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the earth, and thereby obtain the limit of detection. for suffi— 
ciently great distances, we may neglect Noe and o in Eqs. 


(52) and (55), so that 


Ey yh _ “2 2 [ar r' Ja. pir» (75) 


and the integration is over only the regions in the atmosphere near 
the earth's surface where gamma rays and x rays are deposited. 

Substituting Eqs. (4), (15) and.(16) into Eq. (75) and making 
the delta-function approximation for f(t), we find 


y= = { dr' plan (76) 


_ epe ¥ 
Ey gt 7) oE 


where Jy oi is defined by Eqs. (69) and (70). Since the integra— 
2 


tions extend over a region small compared to r' and since the 


I 
integral | = over the important region is large, 


~ «SBE XY 
E, tT) ® oF Jap’? 2 . (77) 


Finally, it is clear that J a. % (r)~ 1, so that 
? 


eBe Y¥ 
sz (78) 


8,9 2E 
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For gamma rays, 


2 
~ 6x 10 
Boy * 2 


esu , (79) 


where the fraction” of the total yield Yip emitted as gamma rays 


is taken to be 10> . For x rays, 


esu , (80) 


where half the total yield” Yop is assumed to be emitted as x rays. 

To estimate the maximum detection range of this signal, we 
shall assume that the only interfering noise is due to the cosmic 
noise background. Thus, it will be assumed that the detectors can 
be effectively shielded from man-made disturbances. 


The mean cosmic noise power is 


22 


= 4nkT/% =1.7x 10. oT w/m* ; (8 1) 


PN 
where T is the effective temperature in degrees Kelvin of the 
cosmic noise at the wavelength \ in meters and o is the detector 
bandwidth in cps. Since the signal to be detected has frequencies from 
one to ten megacycles (See Eqs. (66), (67), (69), (70) and (7 1)), the 
detector bandwidth should be adjusted to o ~ 10! cps at \~ 10° m. 
Under these circumstances the effective temperature of the cosmic 


noise backgroundis Tw~ i0° °K, Hence, the interfering noise power 
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is 


13 2 . (32) 


P._ = 1.7x10 ~“w/m 


N 


For reliable detection—using coincidence techniques—we require 
the peak signal power P to exceed the noise power by a factor of 


ten, that is, 


Cc 2 
— = . 8 
P. ae E. 10 Py (83) 
Or 
—10 —5 
E, * 8 x 10 esu or 2.4x 10 ° v/m. (84) 


Thus, from Eqs. (79) and (80), the limiting detection range for 


gamma rays 1s 


re, @ 0-8 x 10° V Yuen 3 (85) 


and for x rays is 


7 
rv 1-5 x 10 m Yigg (86) 
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VI. DISCUSSION 


An interesting consequence of Eqs. (72) to (74) and Eq. (78) 
to (80) is that for near explosions, the gamma-ray generated electro— 
magnetic signal is more intense than the x-ray generated signal, 
while for distant explosions the opposite is true. To understand this 
result we note that for distant explosions the intensity of the radiated 
field is proportional to the yield of the explosion appearing as gamma 
rays or x rays—thereby, as would be expected, the x-ray generated 
Signal should be larger. As the burst point comes nearer to the 
earth's atmosphere, the intensity of the signal begins to level off 
and becomes independent of yield. For xrays, this leveling off occurs 
at a very much greater distance than for gamma rays—thus the 


smaller intensity of the x-ray generated signal. 


A final point—in concluding that the detection range of the 
x-ray generated signal was given by Eq. (86) and exceeded the gamma— 
ray generated signal, it was assumed that the x-ray yield of the 
explosion was about one—half the total yield. Actually, however, 
this yield might be considerably less—perhaps by a factor of 107 to 
10° , if shielding were employed’. In this case the limiting detection 
range would be determined by the gamma-ray generated signal. 
Shielding of the gamma-ray signal might also be possible, thereby 


degrading the detection range even further. 
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APPENDIX 


In this appendix we derive the formulae for the primary elec— 
tron current and secondary electron density as functions of space and 
time. 

The rate of production of primary electrons at a point r', at 
time t', with initial motion in direction © , is given by 

r' 


A pits REQ) = g(r') f(t -—) h @) , (A—1) 


where the angular distribution h(@2) is so normalized that 
| aa h(Q) = 1, (A—2) 


[a2 a nia) =e, (A-3) 


and 2 is the cosine of the angle between the direction of motion of 
the electrons and the initial gamma-ray or x-ray direction. For 
gamma-ray produced Compton electrons, «* 1. For photoelectrons 


mt 


from xrays, € = v/e ~ 0.05to 0.i for 1ltoa few kev electrons. 
From symmetry, the other components of &. integrate to zero. 
The primary electron current and density at a position r 


and time t are obtained by adding the contributions from those 


electrons made at earlier times and appropriate positions such that their 
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motion brings them to r attime t. That is, 


r' 
(r,t) = ~e Jar' [att [ae g(e') f(t! —)n(a) vit, t':2) 


sf nor - fae" vie e2 | , (A—4) 
i! 
where v (t,t';Q) is the velocity at time t of an electron made at t' 
with initial direction 2. The density of primary electrons is given 
by the same expression without the factor ~—ev. In the absence of a 
magnetic field or any slowing down mechanism, the velocity of the 


- electrons is constant, so that 


v(t,t59) =v 2 , (A-5) 
t 
| dt" ve" t5@) = v a(t—t). (A-6) 


it 
Since the electrons have a finite lifetime given by their range R divided 
by their initial speed, the integral over t' in Eq. (A—4) extends only 
over the finite range from t-— Riv. to t. We shall assume that 

the distances the electrons travel are small compared to r, which is 
the case for currents generated in the atmosphere by high—altitude 
nuclear explosions. Thus, the direction from the source to r' can be 


taken to be the same as that to r and the current simplifies to 


t J 
i@,t)s—ev g(r) [ar' | ate" - 2) [ang nig) 
tRiv. 


e p[ rat voit 1, . (A-7) 
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With t-—t' =T', t—r/c =T, we have 


Riv. vm 
j(r,t) =— ev g(r) J dr [de h(Q) Q £ (1- t+ 2). (A-8) 


For Compton electrons h(@) is essentially a delta function in 


the initial direction. Thus the current is radial and 


Riv v 
j (2) t) =— ev g(t) | Car" f [ -— (1 — = ) ; (A~—9) 
O 
Riv. Vv 
Mri t)= g(r) | dr’ | tT—T'(1 — — |. (A—10) 


For photoelectrons, we expand the function f in powers of 


vor R 
29 <—Q . Thus 
c r c rx 


vit vot , 
— ri _—— _ wt — _ wl _ 
f{(r— fT! + 2) f(r—T') + = 2 f(r—m) +... (A—11) 


The integration over 2 in (A-8) converts Eq. (A—11) into an expansion 
ve 


in terms of which is bounded by ae and which is small compared 
to one at all altitudes of importance. Thus for photoelectrons, 


Riv 
jf, t) #— eev gr) J Pdr f (r=), (A—12) 


Riv, 
nts 1) * ete) | de! f(r— 1). (A—13) 


We note that although different approximations have been used, 


we can write the results in one form for both Compton electrons and 
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photoelectrons 
Riv. vo 
j.(, t) we — eev g(r) | a r(1— — ) ; (A—14) 
O 
nit t)=- j,/ee Vo? (A-15) 
4 “o 
where « #1 for Compton electrons, and « = 5. and the term 
V 
O 


= << 1 in the argument of f should be neglected for photoelectrons. 
In the presence of a magnetic field, the electrons spiral 
about the field so that there are transverse as well as radial currents. 
Moreover, the radial velocity of the electrons is not constant with 
time but varies. Nonetheless, the earlier approximations are 


clearly still valid. Thus, 


Riv 


J (r,t) ~~ CC g(r) J o div (t') tr clo + a ; (A—16) 
Riv | 

Nori tt) © g(t) | Pane] ras + an (A~17) 
0 


where v(t') is the velocity of an electron at time r' after it is made 
and x(t') is its displacement in the radial direction (direction of the 


initial gamma ray or x ray) from its initial position. As before, 


PPE ie pe alee a ee i 
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« 1 for Compton electrons, and we drop x(1)/c in the argument of 


{ for photoelectrons. 


If only a static magnetic fieldis present, the motion of the 


recoil electrons is very simple. The component of velocity along the 


—~ - magnetic field is unchanged and the motion in the plane perpendicular 


to the field is circular with angular frequency 


where 


In terms of spherical coordinates, 


x(T) 


2 _ 2. SinwT 
“| £08 87+ sin @ ——— ; 


2 
v7) Vo | cos 6+ sine cos or , 


vr) =-V sin 9@cos 6(1—coswT) , 


(tT) = V5 sin §9sinwT, 


“> 


(A—18) 


(A—19) 


(A—20) 


(A—2 1) 


(A—22) 


(A--23) 


where @ is the angle between the radial direction and the mag— 


netic field. When these expressions are substituted into Eqs. (A—16) 


and (A—17), we obtain the results given in Section II, Eq. (13) to (16). 


We have neglected several effects in deriving these expressions. 


First, we have assumed that the electrons have a constant speed, and 


? 
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F @ then abruptly stop after an appropriate time, R/ Yot Actually, of 
course, they slow down as they lose energy by ionizing the air. Since 
the rate of energy loss, F(¥), is known for electrons, it can, of course, 


be included. We then have to solve simultaneously the equations of motion 


dE 2dav 

q; = me a =— F(y) , (A—24) 
~ evxB 

d - = F 

gq, (Ymy) = — C — W) Vv: (A~25) 


To do this, we note that a straightforward integration of Eq. (A—24) 
gives 
Y 
dv ._itT (A—26) 


F(y) 2 
O mc - 


where Y=¥ (r=0). Solving this equation for v as a function of tT we 


obtain the electron speed 
1/2 


| v(t} =c [a : | . (A—27) 
| ¥ (7) 


Further, letting v = vu, where ul is the unit vector in the direction 


of the velocity, we find 


du A A B 
ar -_— vme UXBr-wlt)uxy, (A—28) 
or 
du &xB 
wo SR? (A~29} 
where 
v 
= t 1 = eB dt’ = o dv 
@(r) = | dat wr) = —— | Gh) eBe | Fi) (A-30) 
O 7) 
S Thus the electron velocities are the same as found before, with v 


—i_ 


replaced by v(t) =c a ——> 
y (1) 


1/2 
| » and wt replaced by @ (7). 


: 
| 
| 
i 
| 
| 
| 


—34— 


The net effect is to reduce the velocity to zero during the slowing down ¢ 
and to increase the Larmor frequency » as v decreases to one. 

Of greater concern is the question—what magnetic field should 
be used in determining the electron motion? Under conditions where 
the currents generate only weak fields, we need use the ambient mag— 
netic field only. However, under conditions where the explosions are 
not too distant, the currents can give rise to intense crossed E and 
B fields which can significantly affect the motion of the electrons. 

For x-ray photoelectrons, however, the fields even under these condi— 
tions are sufficiently small so that only the ambient magnetic field need 
be used. For Compton—recoil electrons, the fields are large so that 
these self—generated fields must be considered. To estimate the 

effect of these Compton—current generated fields, we note first that 

the electric field produced by the electrons is transverse and in such 

a direction as to oppose the transverse motion of the electrons. The 
generated magnetic field, however, is in sucha direction as to add 

to the ambient geomagnetic field and increase the transverse accelera— 
tion. The intensities of the generated electric and magnetic fields, 
which are equal in magnitude, depend upon the amount of transverse 
current and are time dependent. It is clearly impractical to try to 
solve rigorously the self-consistent problem for the electron motions. 


What we can do first is to observe that the total force on an electron is 
vo 
Fe-e/ E+ xB +B,)| (A-3 1) 


where E and B are the fields due to the currents and Bo is the 


ambient field. Since |E| = |B | , the net transverse force is 
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vV vV 
proportional to E (i — =) plus the transverse component of — xB. 


O 

If the radial electron velocity, Vie» is approximately light speed, 
the effect of the induced fields is greatly reduced. For Compton 
electrons, then, the initial turning is due primarily to the ambient 
field. 

We can next approximately solve the self-consistent problem 
for Compton electrons by taking a simple special case. We let the 
gamma-ray source vary as a delta function in time, and take the 


initial direction of the electron motion at right angles to the ambient 


field. Then the equations of motion are 


a7 (yx) = . ¥- as E(t)y, (A—32) 
F eB e 2 
aq, fy) = Mme x —-— — E(t) (1— =), (A-33) 


where x is the direction of initial motion and y is the transverse 
direction. For not-too—distant explosions, for which the fields are 


highest, we find from Eqs. (61), (A—16) and (A—17) that 


y(t) (A—34) 


It is straightforward to show that the maximum field strength occurs 


at t = 0 and has the value 


YB 1 
E(0) =B. ae(l—s) yk (A—35) 
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If we neglect the generated field in Eqs. (A-32) and (A—33), and 
consider the motion only in the ambient geomagnetic field, we find 


the peak value of the electric field again occurs at tT = 0 and is 


given by 
VR 
E(0) =-Bo vac(1—8) . (A—36) 


Thus the effect of the generated fields is to decrease the peak field 


by the factor 


VR —1 
fy + sae | 20.7, {(A—35) 


where parameter values corresponding to 1 mev electrons have been 
used. 

This example indicates that as long as the peak radiated field 
strength occurs near the beginning of the current pulse, the effect 
of the generated fields on the electron motion does not significantly 
change the order of magnitude of the estimate for the peak radiated 
field. The generated fields do affect the time variation of the 
radiated field. But since the range of detection depends principally 
on the peak field strength, the effect of the generated fields on the 


electron motion may be disregarded. 
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